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ABSTRACT 


This  theory  permits  the  calculation  of  pressure 
distributions  over  a  thin  airfoil  in  steady,  inviscid, 
incompressible  flow;  or  given  a  desired  chordwise  pressure 
distribution  the  camber  line  of  an  ideal  wing  can  be 
determined.  The  method  treats  the  circulation  about  the 
wing  -as  a  continuous  vortex  sheet  of  variable  strength 
covering  the  wing  planform  and  trailing  downstream  to  infinity. 
The  results  of  the  present  numerical  solution  for  the  pressure 
distribution  solution  are  not  satisfactory.  However  the 
solution  for  the  camber  line  of  an  ideal  wing  are  in  reasonable 
agreement  with  published  two  dimensional  results. 
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I.  INTRODUCTION 


This  analysis  is  a  continuation  of  two  dimensional 
thin  airfoil  theory  ["Ref.  l]  and  lifting  line  theory  [Ref.  2] 
into  a  three  dimensional  lifting  surface  theory,  using  the 
concept  of  a  continuous  vortex  sheet  over  the  wing  which 
trails  off  to  infinity  in  the  downstream  direction.  The 
continuous  vortex  sheet  differentiates  this  method  from 
the  vortex  lattice  kernel  methods  of  Multhopp  [Ref.  3]  , 
Falkner  [Ref.  4]  ,  Weissinger  [Ref.  5]  ,  and  Watkins  et.  al. 

[Ref.  6]  and  double- lattice  finite  element  methods  of 
Lopez  and  Shen  [Ref.  ?]  and  Giesing  et.  al.  [Ref.  8  ]  . 

The  present  method  is  restricted  to  steady,  inviscid, 
incompressible  flow  and  thin  airfoils  which  are  symmetrical 
in  semi span. 

The  circulation  about  the  wing  is  treated  as  generated 
by  a  continuous  vortex  sheet  of  variable  strength  covering 
the  wing  planform  and  trailing  downstream  to  infinity.  The 
vortex  sheet  strength  is  set  by  the  Kutta  condition  at  the 
trailing  edge  and  requiring  no  flow  through  the  wing  at 
specified  control  points. 

The  flow  equations  for  the  wing  are  resolved  into  down- 
wash  angles,  with  assumed  series  solutions  for  the  circulation 
functions .  From  the  known  planform  geometry  and  the  required 
solution  at  the  control  points,  no  flow  through  the  wing,  the 
constants  of  the  circulation  functions  series  are  determined. 
The  equations  are  put  into  matrix  form,  the  matrices  being 
formed  by  numerical  summation  over  the  planform,  and  solved 
by  matrix  algebra.  Once  the  constants  of  the  circulation 
functions  are  known  all  of  the  airfoil  characteristics  can 
be  calculated. 

This  analysis  generates  the  spanwise  and  chordwise 
pressure  distributions  over  the  wing  which  are  necessary  in 
designing  the  wing  structure.  The  solution  is  also  done  in 
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reverse,  that  is  the  pressure  distribution  over  the  wing  is 
specified  and  the  shape  of  the  mean  camber  line  of  an  ideal 
airfoil  is  determined. 

This  analysis  uses  the  consepts  of  additional  and 
reference  lift  which  is  of  academic  interest  in  under¬ 
standing  airfoil  theory. 
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II .  General  Problem  Development 


A.  THE  WING  VORTEX  SYSTEM 

The  present  analysis  uses  the  theory  of  thin  wings  to 
describe  the  characteristics  of  lifting  surfaces  in  steady, 
inviscid,  incompressible  flow.  The  most  common  planforms 
are  those  which  involve  straight  leading  and  trailing  edges, 
as  illustrated  in  Fig.  1,  but  the  present  analysis  applies 
also  to  wings  having  curved  leading  and  trailing  edges .  We 
consider  only  wings  that  are  symmetrical  with  respect  to  the 
x  axis.  We  also  restrict  the  present  analysis  to  lift 
distributions  which  are  symmetrical  with  respect  to  the  x 
axis . 

The  wing  and  the  trailing  vortex  system  associated  with  it 
are  treated  as  a  continuous  vortex  sheet  of  variable  strength 
y  .  The  sheet  strength  'S  is  regarded  as  a  vector  in  the 
x,y  plane.  The  wing  vortex  sheet  is  treated  as  lying  in  or 
very  close  to  the  x,y  plane  and  &  is  everywhere  tangent 
to  the  vortex  sheet. 

The  sheet  strength  <T  is  related  to  the  velocities  Vu 
and  Vl  just  above  and  just  below  the  vortex  sheet  at  a 
given  point  x,y  and  to  the  remote  velocity  \/oo  ,  as  indicated 


in 

the  vector  diagram,  Fig.  2. 

From  the  vector  diagram  Fig. 

2. 

Vu  =  \7eo-  i  Y*  £ 

(1.1) 

Vl  =  V«  -  J  Y  x  I 

(1.2) 

It 

then  follows  that 

( 1/u  Vl  )  =  Vo o 

(1.3) 

(l/u  ~Vu)  = 

(1.4) 
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or  conversely 


y=  k*  ( Vu-  Vu) 


(1.5) 


In  Fig.  3,  consider  the  circulation  d  T  around  contour 
AOBB  0  A' A.  Line  segment  AOB  lies  just  above  the  vortex 
sheet,  and  segment  B/Q/A/  lies  just  below  the  vortex  sheet. 


AOs  =  A  O  S  =  d 


(1.6) 


Then 


d  r  c  ^  V'  dr  -  V  9  J  r  V*Jr  +  J  l /*  d  r  +  i  l/*J 

C  A  4  B  A' 

d  T  -  l/u*Js  +  0  -  Vt-’ds  +  0=  (  Vw  -  Vt_  )  •  d  s 

dr  ^  (f*f)  •  js 


r 


(1.7) 


It  is  assumed,  which  will  subsequently  be  confirmed,  that 
the  dT  in  equation  (1.7  )  is  an  exact  differential.  There¬ 
fore,  there  exists  a  function 

r--  ru.a) 

such  that  for  any  small  changes  in  the  coordinates 


df-  (w-) +  (tv) 


(1.8) 


Translating  this  into  vector  terms  gives  : 


«jr=(t^  +  T-|ij-(rdx*jdS)=vr-Js  a.?) 


bx  ■  -  by 

Comparing  equations  (1.7)  and  (1.9)  shows 

V  r  =  V  < 

or  conversely  that 

V  =  k  x  vr 


(1  .10) 
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(1.11) 


Fig.  1  Wing  Coordinates 
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y 


i'ig.  2  Velocity  Vectors 


rig.  3  Vortex  Sheet  Element 
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These  results  confirm  that  is  an  exact  differential. 

Translating  equation  (1.11)  into  Cartesian  components 
gives 


(1  .12) 

(1.13) 


The  vector  function  V  defined  by  equation  (1.11)  or  by 
scalar  equations  (1.12)  and  (1.13)  is  non-divergent . 


7 9  V  - 


(1 .14) 


B.  LIFT 

Newton’s  Law  applied  to  a  fluid  element  of  fixed  mass 

pAXAyA Z.  is 


F 


^  Ax  Ay  Az. 


oy 

£>t 


(2.1) 


It  is  assumed  that  the  fluid  is  non  viscous;  that  is,  shearing 
stresses  are  absent.  The  resultant  pressure  force  of  the 
fluid  element  arising  from  the  static  pressure  variation  in 
the  fluid,  is  the  negative  of  the  pressure  gradient  multiplied 
by  the  increment  volume. 


resultant  =  (  -  grad  p  ) 


The  pressure  force  plus  the  weight  is  the  external  force, 
therefore  the  equation  of  equilibrium  becomes 

f  3  -  SrQ'cf  P  -  ?  (2*2) 

which  is  Euler's  equation. 

If  the  fluid  is  assumed  incompressible  and  equation  (2.2) 
is  integrated  dropping  the  unsteady  term  it  becomes 
Bernoulli ' s  equation 

\/z 

?4r  ’P- 


constant 
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(2.3) 


The  pressure  difference  across  the  wing  at  an  arbitrary 
point  can  be  found  from  Bernoulli's  equation. 


/  PL-  Pu 

\  -  , 

/'API 

)Z_ 

.  /- 

Vu 

l  <£co 

/  ~  1 

^  ‘ 

^  |/«  J 

l/oo 

 (  Vu  + 

Vl.)  • 

(  Vo  -  Vu ) 

2  i/L 

•  ( 

y* 

t) 

\4> 

2. 

vz 

 Z'i  • 

(  V 

X  k.  ) 

r 

2 

V.. 

l/«. 

V- 

(2-4) 

Equation  (2.4)  ultimately  fixes  the  lift  distribution  over  the 
wing  surface  once  the  circulation  function  P  has  been  found. 

In  the  trailing  vortex  region  behind  the  wing,  there  is 
no  wing  surface,  and  hence  no  pressure  difference  across  the 
vortex  sheet.  According  to  the  Kutta  condition,  the  pressure 
difference  drops  to  zero  not  only  behind  the  wing  but  also  all 
along  the  trailing  edge.  This  is  expressed  as  follows: 


(2.5) 


This  result  shows  that  the  vortex  lines  behind  the  wing  and 
along  the  trailing  edge  are  all  parallel  to  the  x  axis,  with 
no  spanwise  components.  Spanwise  vorticity  "jf  can  occur 

only  along  the  wing  surface. 

/ 

In  Fig.  4,  next  page,  consider  the  incremental  lift  d  L 
exerted  by  the  doubly  shaded  portion  of  the  wing,  that  lies 
between  y  and  (y+dy)  and  between  x  and  x.  Thus: 

d  L  =  i  v  V*,  J"  (  )  dx  d'd 

Xu  ^ 

O 

d  L  =  £  14*  r  dy  (2.6) 


This  result  is  a  form  of  the  Kutta-Joukowski  Law  which  states 
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chat  the  force  experienced  by  a  body  in  a  uniform  stream  is 
equal  to  the  product  of  the  fluid  density,  stream  velocity, 
and  circulation  and  has  a  direction  perpendicular  to  the 
stream  velocity.  Here  P  represents  the  circulation  about  a 
contour  C  in  the  Plane  y  =  constant.  Contour  C  passes  in 
front  of  the  wing  and  pierces  the  wing  surface  at  point  x,y. 

If  we  now  extend  the  upper  limit  of  integration  to  the 
trailing  x^,  holding  y  constant,  we  obtain: 
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X  X  T 

r  ( *>  y)  n;  (y) 

ciL'  ,  JL 

dL  -  e  (/» r; dy  (2.7) 

thus  fr  (y)  fixes  the  spanwise  lift  distribution  over  the  wing 
In  order  to  reduce  equation  (2.7)  to  a  convenient  dimension 
less  form,  it  is  useful  to  introduce  the  following  auxiliary 
notation . 


s  = 

total  wing 

2 

area,  ft 

(2.8) 

2s  = 

b  =  wing 

span,  ft 

(2.9) 

c  = 

=  mean 

chord,  ft 

(2.10) 

Ak  = 

II 

”ji|^ 

II 

aspect  ratio 

(2.11) 

n  = 

M 

—  =  spanwise  coordinate 

(2.12) 

ii 

dy 

s 

(2.13) 

d  1— 

section  lift  coefficient 

(2.14) 

o^«.c 

11 

J 

<U 

L 

-S 

wing  lift  coefficient 

(2.15) 

By  utilizing  the  above  nomemclature,  it  is  easy  to  reduce  the 
spanwise  lift  distribution  as  given  in  equation  (2.7)  to 
either  of  the  following  two  equivalent  forms. 

(2.16) 

To  find  the  total  lift  we  integrate  over  the  span,  For 
symmetrical  lift  distribution  it  is  only  necessary  to 
integrate  over  the  semi-span  and  multiply  by  two.  Thus: 
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(2.17) 


O  O 

C.  INDUCED  VELOCITY 

The  velocity  cWP'  induced  at  an  arbitrary  point  P, 

coordinates  x  ,y  on  the  wing  surface  by  an  element  of  the 
P  P 

vortex  sheet  which  lies  at  the  point  whose  coordinates  are 
x,y  is  found  as  follows.  The  vortex  sheet  element  has  length 
ds,  width  dn  and  strength  Y  .  The  circulation  strength  of 
the  vortex  filament  involved  is  =  Ydn. 

The  Biot  Savart  Law  for  this  case  may  be  written  in  the 
form : 


dvv  p 


IT x  Y  dn  ds 
4  tt  r  3 


(3.1) 


where 


T  =  i  (x-x?)  +  -j  (  4  -  HpW  k  Cz-"zp) 


(3.2) 


dn  ds  =  dS  =  element  of  the  wing  area 
y-  c  v*  k  Yz.  -  k  *  V  T 


(3.3) 


Neglecting  the  terms  in  equation  (3.2)  and  k.  V"z.  in 

equation  (3.3). 


also 


r  4 


i 


i 


k. 


(x-xp)  I'q-'dr) 

-  —  vr 

&  'j 


o 

o 


-  k  [(*-x^(^r)  +  (  1 


(3.4) 


(3.5) 
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Y 


A 


Fig.  5  Induced  velocity  in  the  x-y  plane 
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t  t 


Upon  substituting  expressions  (3.2),  (3.3),  and  (3.5)  into 
(3.1),  and  integrating  over  the  entire  wing  surface  S 
(  but  not  yet  over  the  vortex  sheet  trailing  behind  the  wing) 
the  downwash  due  to  the  wing  vortex  is  obtained. 


W  P  = 


!e 

TF 


'  £(*-*»)  TV  ^  3 

s  [(x-*,)*-.  ]iA 


(3.6) 


The  incremental  downwash  velocity  dvvp  induced  at  point 
P  due  to  the  semi-infinite  trailing  vortex  filament  of  strength 
cl  Hr  which  leaves  the  wing  at  the  point  whose  coordinates  are 
xT  and  y  as  shown  in  Fig.  5  is  by  the  Biot  Savart  Laws 


y 


I 


dv;  =  jl.  ( 

4TT  '  y-y? 


(X-Xp) 


*  k 


(X-X&) 


-$> 


(_L_\  ,  _ 

[(x-x.OVv'dpf  ] 


5^ 

(3.7) 

the 

wing  gives 

• 

1 

H 

(3.8) 

The  total  induced  velocity  at  point  P  is  the  sum 


WP  =  M/p  -+■  Wp  (3.9) 

where  Wp  and  vvp  are  evaluated  respectively  from  equations 
(3.6)  and  (3.7).  The  first  of  these  integrals  represents  the 
velocity  induced  by  the  vorticity  over-  the  actual  wing  surface, 
and  the  second  integral  represents  the  velocity  induced  by 
the  vorticity  trailing  behind  the  wing. 

The  downwash  integrals  are  now  non-dimensionalized  and 
the  Cartesian  coordinates  x  and  y  are  transformed  into 
dimensionless  relative  coordinates  f  and  £  defined  as 
follows  s 
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=  relative  chordwise  coordinate 


(3.10) 


X-Xt_ 

C- 

d  n  . 

j-  =  relative  spanwise  coordinate 

The  following  auxiliary  definitions  are  used. 

c  =  -j—  =  average  chord 

/q  Zs. 

™  =  — —  =  aspect  ratio 

(?)  =  — r—  =  leading  edge  function 

2!r(7)  =  =  trailing  edge  function 

1  c 

=  C  (7)  =  relative  chord  function 

C<.?)  =  Xt(T)  -2.J?) 


(3.11) 

(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 


Note  that  if  the  parameter  is  specified  and  if  any  two 
of  the  three  functions  involved  in  equation  (3.17)  are 
specified,  this  information  completely  defines  the  wing 
planform. 

The  circulation  function  T  is  non-dimensionalized  as 
follows.  The  number  four  is  inserted  into  these  definitions, 
somewhat  arbitrarily,  because  this  reduces  the  final  downwash 
integrals  to  a  particularly  convenient  numerical  form. 

=  dimensionless  circulation 
function  along  the  wing 
surface  (3.18) 

=  dimensionless  circulation 
function  at  and  behind  the 
trailing  edge  (3.19) 


r *(}>•?)  - 


r.d)  = 


Tt 

”4  Vco-s 
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The  induced  velocities  Wp  and  vvr  are  non-dimensionalized 
to  corresponding  downwash  angles  c<^  and  «:P"  as  shown  below. 
Those  angles  are  defined  as  positive  if  down  from  free  stream 
velocity . 


/ 

0 **  p 


ft 

(Xp 


Wp  « la. 


downwash  angle  point  P  induced 
by  vorticity  distribution  over  the 
wing  (3.20) 


— ft  . 

Wo  «  fe. 
V'oo 


downwash  angle  at  point  P  induced 

by  vorticity  distribution  behind 

the  wing  (3.21) 


p  =  total  downwash  angle  induced  at  point 
P  by  the  complete  vortex  sheet  system 
of  the  wing  (3.22) 


Some  useful  geometrical  relationships  are  : 


(X-Xp)  =  s,  JL[(Xl-Q)  -(^up  +Cp?p)  ~]  (3.23) 

(oj-^p')  -  s  (  ^  -  *?p)  (3.24) 

Also  at  any  point  on  the  wing 

r(x^)  ^  "M/oos  r* ( ?  >  y)  (3.25) 

Differentiating  equation  (3.25) 


Differentiating  equations  (3.24)  and  (3.25) 
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=  s,  -^r  [_C  A  l  +  (X!l  +  C  f  )  d  *X  "  °  3 


c{x  -  -  ^ 
dy  =  s.  d  >1 


(3.27) 

(3.28) 


The  prime  marks  in  equation  (3.27)  denotes  differentiation 
with  respect  to  y.  Substituting  equations  (3.27)  and  (3.28) 
into  equation  (3.26)  gives 


=  (t£)  4  +  C  1 


c>  *1 


ft 


(3.29) 


The  coefficients  of  dz  and  dy  must  be  separately  equal  on 
both  sides  of  equation  (3.29)  since  this  relation  must  be 
satisfied  for  arbritrary  values  of  dz  and  dy.  Hence 


(£)*  c  = 

) 

(3.30) 

/  ^ 

)  ^ 

(3.31) 

^  r 

Solving  those  two  equations  for 

V  P 

and  - —  gives 
^  & 

* 

M 

(3.32) 

(&)- 

(3.33) 

The  element  of  wing  area  in  relative  coordinates 

may  be 

expressed  in  the  following  manner. 

dS  = 

>*-s 

dS--  Zsz  (z\  )(?  )  = 

(3.34) 
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dS 


Fig.  6  Element  of  Wing  Area 


The  various  quantities  of  equations  (3.6)  and  (3.8)  are 
transformed  into  non-dimensional  notation. 


f)1 1 


ir 


=  t'C*’ lg-[(XsC[)-(Xa-+C,U)]  {{y 


d  l 


-  W-'M  C(x^-)  dfc!*  (3.36) 


(xT  -X»)  =  s.  ^  [r-  -  *  Cr  !p) 


(3.37) 
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Study  of  the  algebraic  structure  of  these  relations 
shows  that  the  final  downwash  integrals  can  be  reduced  to  a 
relatively  concise  form  by  adopting  the  following  auxiliary 
nomenclature.  Let 


[~z[(Xl+Q.[)-(^p<f!p)1Z+  (l-tlr)*  J*z 


F*  = 


-  0- 


=>5 


(3.39) 


(3.40) 


The  downwash  angle  integrals  equations  (3.6)  and  (3.8)  can 
now  be  expressed  in  the  following  form. 


(3.41) 


(3.42) 


D.  THE  CONDITION  OF  NO  FLOW  THROUGH  THE  WING 

The  shape  of  the  mean  wing  surface  near  the  x,y  plane 
may  be  expressed  by  the  function 


z  =  z(x,y) 


(4.1) 


A  unit  vector  nV  normal  to  the  above  surface  at  the  arbitrary 
point  P  on  that  surface,  in  terms  of  its  Cartesian  components, 
is  given  by  the  expression 
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_  -(jit)?  i.  -  (,4h*  i  *l!lk 

n"=  U¥)^i^TTX^--  (4. 

Since  and  (-5^-)^  are  negligible  in  comparison  with' 

unity,  equation  (4.2)  is  adequately  approximated  by  the 
linearized  version 


ho  - 


=  - 


'b  Z 


hu)  )p  \ 


a)  ^ 


(4.3) 


The  resultant  velocity  vector  at  point  P  may  be  closely 
approximated  by  the  expression 

Vp  -  T  \4>  -*•  k.  l/c*  SIM  X  -  )z  \AJp  (4.4) 

Equation  (4.4)  implies  that  the  induced  velocity  w/p  is  taken 
as  positive  if  down. 

Since  x  is  a  small  angle  equation  (4.4)  may  be  linearized 
to 

Vr  rz  i  V<*>  +  k  (4.5) 


The  condition  of  no  flow  through  the  wing  now  requires  that 


.  n  -  o 


(4.6) 


thus 


=  o  (4.7) 

(l4>*-WP)  =  0  (4.8) 


hence 
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tlz.\  - 

^  'bx  - 


<K  - 


Wp 

V5 * 


(4.9) 


but 


*<< 


V 

t*  CC.p 


therefore 


(4.10) 


The  overall  wing  equation  is  now  found  by  substituting 
the  integrals  of  equations  (3.41)  and  (3.42)  into  equation 
(4.11).  The  result  is  the  basic  wing  equation. 


The  influence  functions  F^  ,  F^,  and  F^  are  as  defined  by 
equations  (3.38),  (3.39)  and  (3.40).  These  three  functions 
depend  only  on  the  planform  and  aspect  ratio  of  the  wing. 
Theoretically,  equation  (4.12)  must  be  satisfied  at  all  points 
in  the  wing  surface  for  an  exact  solution.  However,  an 
adequate  approximate  solution  can  frequently  be  obtained 
by  satisfying  equation  (4.12)  at  a  sufficient  number  of 
discrete  control  points  suitably  distributed  over  the  wing 
surface . 

E.  REFERENCE  AND  ADDITIONAL  LIFT  DISTRIBUTIONS 

Let  the  wing  slope  function  be  known  at  every 

point  ,  tjp  of  the  wing.  Also  let  the  entire  leading 
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edge  of  the  wing  lie  exactly  in  the  x-y  plane  and  the 
trailing  edge  lie  in  the  x-y  plane  at  midspan.  In  general, 
however,  z  is  not  necessarily  zero  at  other  points  along 
the  trailing  edge.  Let  17*  (t,?)  be  a  function  which 
satisfies  equation  (4.12)  when  has  a  particular  value  . 
The  subscript  r  stands  for  reference  condition.  Thus  at  the 
reference  condition  equation  (4.13)  becomes 


The  general  solution  or  eQuation  (4.13)  can  now  be 

written  as  the  sum  of  the  reference  solution  and  an  additional 
solution  in  the  form 

ni.'?)  =  (5.2) 


To  demonstrate  that  equation  (5.2)  is  valid,  and  to  obtain 
the  equation  that  governs  the  form  of  the  additional  solution 
function  we  proceed  as  follows.  Substitute  equation  (5.2) 
into  equation  (4.12).  Subs tract  equation  (5.1)  from  the 
result.  Finally  divide  through  by  the  factor  (^oc- ccr  ).  The 
result  is 

_y  [f,  (ifv 

-i  ° 


(3  - 


TT 


+ 


r 


j 

-i 


cj  fr 


V 

To- 


d  v\_ 


)  Jf 


(5.3) 


Thus  equation  (5.1)  governs  the  reference  solution  or  the 
particular  solution  while  equation  (5.3)  governs  the 
additional  solution.  The  specific  value  of  in  equation 

(5.1)  is  fixed  by  the  particular  condition  that 
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z  =  0  at  |=  +1,^=0 


(5.4) 


given  that 


z  =  0  at  £  =0  for  all  f  (5.5) 

The  reference  solution  depends  not  only  on  the  wing 
planform  and  aspect  ratio  but  also  on  the  wing  slope  function 
and  the  angle  of  attack  .  However,  the  additional  lift 
function  involves  only  the  wing  planform  and  aspect  ratio. 

It  will  be  shown  later  that  the  wing  lift  curve  slope 
depends  only  on  the  additional  lift  function  n*T  and  is 
independent  of  the  reference  solution  Tp  •  Therefore, 
the  wing  lift  curve  slope  depends  on  wing  planform  and 
aspect  ratio  only,  and  is  independent  of  the  wing  slope 
function. 

To  specify  the  design  of  a  wing  it  is  necessary  to  fix 
two  planform  functions  and  c(l)  »  the  aspect  ratio 

AK  ,  and  the  wing  slope  function  (-57-)  p  at  all  points 

r[P  on  the  wing.  For  a  wing  of  specified  design,  equation 
(5.1)  must  be  solved  for  the  initially  unknown  function 
fV  and  equation  (5.3)  must  be  solved  for  the  initially 

unknown  function  pp  (1  j?)  • 

In  the  case  of  equation  (5.1)  the  above  procedure  can 
also  be  partially  reversed,  that  is,  the  function  pp  {]>*[) 
can  be  specified  arbitrarily,  within  certain  limits,  and  the 
equation  can  then  be  solved  for  as  the  unknown.  This 

reversed  solution  procedure  had  two  advantages.  Firstly,  it 
means  that  the  lift  distribution  over  the  wing  in  the 
reference  condition  can  be  stipulated  independently  and  that 
the  wing  slope  function  can  always  be  found  such  as  to  yield 
the  desired  reference  lift  distribution.  Secondly,  it 
happens  that  equation  (5.1)  is  much  easier  to  solve  when  the 
wing  slope  function  is  the  unknown  than  when  the  reference 
lift  distribution  is  the  unknown  because,  in  the  latter  case, 
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the  unknown  function  is  under  the  integral  sign.  If 


is  the  unknown,  equation  (5.1)  can  be  solved  by  direct 
integration;  usually  this  must  be  performed  numerically. 
Repeated  integrations  of  equation  (5.1)  are  required,  a 
separate  numerical  integration  being  needed  for  each  discrete 
point  at  which  the  wing  slope  function  is  being  evaluated. 

In  order  to  represent  the  wing  slope  function  adequately 
it  is  necessary  to  choose  a  sufficient  number  of  discrete 
points  suitably  distributed  over  the  wing. 

Unfortunately,  the  above  reversed  solution  procedure 
cannot  be  applied  to  the  additional  lift  distribution  ToT 
as  governed  by  equation  (5.3).  Of  course,  the  equation  can 
still  be  solved,  at  least  approximately. 

The  choice  of  the  reference  angle  of  attack  <xr  which 
fixes  the  corresponding  reference  solution  Vr  is  to  a  certain 
extent  arbitrary,  The  most  common  choice  is  to  set  c<v  equal 
to  the  angle  of  zero  lift  of  the  wing,  that  is  to  the  angle 
which  yields  a  wing  lift  coefficient  of  zero.  For  a 
symmetrical  lift  distribution  the  lift  coefficient  is 


(5.6) 


o 


The  restriction  that  the  lift  coefficient  equal  zero  does 
not  necessarily  mean  that  the  integrand  is  identically  zero 
all  along  the  span,  just  that  the  net  area  inder  the  curve 
sums  to  zero.  Thus  at  =  0,  portions  of  the  span  may  be 
exerting  upward  forces,  while  other  portions  are  exerting 
downward  forces .  A  wing  that  behaves  in  this  way  at  =  0 
may  be  said  to  be  aerodynamically  twisted.  Conversely  a 
wing  that  at  =  0  yields  c^-^-  identically  zero  for  all 
values  of  ^  may  be  said  to  be  aerodynamically  untwisted; 
this  however  does  not  imply  the  absence  of  camber. 

The  reference  condition  for  this  analysis  is  that  which 
corresponds  to  the  angle  of  attack  for  which  the  wing  lift 
coefficient  equals  zero.  This  particular  lift  distribution 
is  also  commonly  labeled  as  the  basic  lift  distribution. 
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F.  ANGULAR  COORDINATES 


A 

?>< 

made 


then 


further  shift  from  the  linear  relative  coordinates 
to  corresponding  angular  coordinates  0  ,  e  ,  is  now 

See  Fig.  7.  The  basic  conversion  relations  are 


2  - 

h  (  i 

-  COS  P  ) 

(6.1) 

f  = 

cos  e 

(6.2) 

oil 

h  sin 

/ 

p 

(6.3) 

r 

ii  _ 

d© 

-  sin 

© 

(6.4) 

n 

)  = 

2  ^  (  \ 

(6.5) 

—  ( 

)  =  ' 

'  —  (  'l 

Si^sl  0  ^  3 

(6.6) 

Fig.  /  Angular  Coordinates 
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The  downwash  integrals  can  now  be  converted  as  follows 


-I  o 


O  TT 


—  iKp  — 


TT/3^ 


TT  O 


TT  TT 


—  - 


T TP^ 


ft  s.m r2  (HI) 

n  ^  J  2  2.  V.  b  ©  J 


6  4  d 


o  o 


TT  IT 


|  j" f°"  +  )  i d?<  d 


O  o 


-«r=w[F>(0H- 

TT 


-  X  p  - 


where 


_L 

tt 


*  (^)  ^ 


(6.7) 


(• 

(6-8) 


Gi  =  sin  6 


Go  =  -  S±  sin 


G3  "  "  F3 


(6.9) 

(6.10) 
(6.11) 


Now  by  refering  to  the  above  relationship  and  to  the  previous 
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equations  which  define  the  functions  ,  and  F^  it  is 

possible  to  summarize  the  final  wing  equations  as  the 
culmination  of  a  series  of  systematic  calculation  steps. 
The  following  scheme  of  auxiliary  notation  is  used. 


l  = 

h  (  1  -  cos  <f>  ) 

(6.12) 

Ip  = 

hi  1  -  COS  ) 

(6.13) 

7  = 

cos 

(6.14) 

= 

cos  ©p 

(6.15) 

(it  = 

iu(eV  C-CM I 

(6.16) 

hip  = 

Xu  (_  ©p)  +  C 

(6.17) 

hi  = 

■i  Lz>-  - 1 L  Lc- (l 

(6.18) 

Depending  on  planform,  the  derivatives 

d  © 

and  i£. 
d  © 

may  be  multi-valued  at  midspan,  however 

for  a 

symmetrical 

planform  they  are  set  to  zero. 

h3  =  2 

sin 

Hip  )  -  i>\-Kp)  hz  ] 

(6.19) 

[Li 

Z  z  ' 

?  3/z 

h4  =  ! 

(ht-  hip  )  ]  -f  ^p) 

[ 

(6.20) 

Sr=  XL-* 

c 

(6.21) 

h5  =  ■ 

(6.22) 

-  7 

,2  7  '/2. 

h6  “  \ 

^  hs 

*  (  7  -  tr)  J 

(6.23) 

G1  “ 

hs 

TS 

(6.24) 

= 

U 

f-'lp)  C 

(6.25) 

2 

h4 
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Os  -  — 


(6.26) 


The  basic  wing  equations  become  the  following 


Reference  lift 


IT  tt 

r 


=  °<r  +  if*.  f  |LG'1  ^ 


■  Te  i  j 


© 


(6.27) 


Additional  lift 


O  - 


c> 


G.  APPROXIMATE  SERIES  SOLUTION 
The  reference  and  additional 
represented  by  double  series  of 


(6.28) 

distributions  can  be  adequately 
the  following  form. 


W  £ 


fr  (  <^)  ©  )  -  Sr  s  H  r  (<3$ ) 


Oro 

ftl  s 


(7.1) 


(7.2) 


For  a  spanwise  symmetrical  lift  distribution,  only  odd  values 
of  s  are  used.  Ref.  1.  There  are  (R-rl )  distinct  functions  Hr 
and  (  ^ )  distinct  sine  functions  in  the  above  series. 

The  total  number  of  initially  undetermined  constants  in  each 
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of  the  above  series  is  therefore 

T  =  (  R  +- 1  )•(  )  (7.3) 

If  there  are  (R+l )  chordwise  stations  and  (  ■S>'gy-  )  s panwise 
stations  over  the  semi-span,  this  establishes  T  control 
points.  By  satisfying  each  governing  equation  at  each  of  the 
T  control  points  we  can  solve  for  each  set  of  T  constants . 
Once  the  constants  Brs  and  Ars  are  found,  all  aspects  of  the 
wing  performance  can  be  readily  calculated  including  lift 
distributions,  slope  of  the  lift  curve,  induced  drag,  etc. 

The  character  of  the  H  functions  is  established  as 
follows 


AP  _  2  ( 

VcoK**) 

(7.4) 

>> r  zLfi^v<*>  fir* \ 

CL  v  > 

(7.5) 

therefore 

AP  /4-fc\  (>>?*)  (  4-/^  ' 

\  2. 

(>r*) 

(7.6) 

'  CL  ~  ^  C_ 

/  Sifvi  9$ 

For  the  additional  lift,  equation 

(7.2) 

is  substituted 

into 

equation  (7.6) 

C  Cr  O 

r  2 

1  S/M  S  © 

M0  H 

(7.7) 

For  the  reference  lift,  equation 

(7.1) 

is  substituted 

into 

equation  (7.6) 

R+i  -s 

(f-)r=0  HBr^ 

z 

SiM  s  © 

>(Ni 

l\ 

(7.8) 

It  is  now  possible  to  make  use  of  the  results  of  thin 


airfoil  theory  and  set  up  the  functions  in  equations  ( 
and  (7.8)  in  an  analogous  form.  Ref.  2.  Thus  for  the 
case  of  r  =  0  we  use 

4  __  2L  _  \  4-  cos  <& 

'b  |  d  ~  S/M  jz5 

and  for  r  ^  0  we  use 


Equations  (7.9)  and  (7.10)  can  be  integrated  to  yield 
For  r  =0 


Ho(^  =  4r  (  +  S/M  J* 'l 

For  r  =  1 

H, 

For  r  >  2 

|_j r  =  JL  f  siNj(r-i)  _  sim  (r-vi)  gS  7 
r  l  r- 1  r+ 1  J 

At  the  trailing  edge,  0  =  it  ,  these  relations  yield 

HoM-  TL 

HAtt)  =  ~ 

Hr-(TT)  -  O 


.7) 

special 

(7.9) 

(7.10) 

(7.11) 

(7.12) 

(7.13) 

(7.14) 

(7.15) 

(7.16) 


By  substituting  these  values  into  equations  (7.1)  and  (7.2) 
the  corresponding  expressions  for  the  trailing  vorticity 
are  obtained. 


(7.18) 


fn-a.  L&)  ~  ( 2-  Aoi  +  S  "J  5© 

S»),^,?... 

An  ideal  wing  is  one  for  which  =  0  at 

the  leading  edge.  This  requires  that  Bos  =  0  for  all  s. 

If  the  planform  function  is  arbitrarily-  prescribed  the  wing 
should  not  be  presumed  to  be  an  ideal  wing  and  all  terms  in 
Bos  should  be  included  in  the  solution.  The  additional 
lift  is  always  characterized  by  the  singularity  (•^2)— 00 
at  the  leading  edge  so  that  the  terms  in  Aos  must  always 
be  retained. 
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Ill  DESCRIPTION  OF  NUMERICAL  PROBLEM  SOLUTION 


A.  CALCULATION  SCHEME 

The  calculation  scheme  in  this  analysis  uses  two  distinct 
meshes.  There  is  a  coarse  mesh  of  control  points,  and  a 
fine  mesh  for  the  purposes  of  accurate  integration.  It  is 
desirable  with  respect  to  both  meshes  that  the  ratio  of  the 
number  of  spanwise  stations  over  the  semispan  to  the  number 
of  chordwise  stations  should  be  roughly  equal  to  the  aspect 
ratio  divided  by  two.  This  gives  roughly  equal  resolution 
in  the  spanwise  and  chordwise  directions. 

The  total  angular  interval,  0  to  n-  in  each  case  is 
subdivided  into  equal  sub-intervals  and  the  calculation  point 
located  at  the  center  of  each  sub-interval.  The  calculation 
points  associated  with  the  fine  mesh  must  not  coincide  with 
any  of  the  control  points  associated  with  the  coarse  mesh 
because  this  would  cause  the  functions  and  F2  to  assume 
an  indeterminate  form.  While  this  indeterminacy  can  be 
resolved  it  is  best  avoided.  A  good  way  to  avoid  this 
difficulty  is  to  make  the  fine  mesh  an  even  submultiple  of 
that  of  the  coarse  mesh.  Then  the  control  points  will 
always  lie  along  intersections  of  the  boundaries  of  the  area 
elements  which  comprise  the  fine  mesh.  See  Fig.  8. 

The  double  integration  over  the  wing  area  reduces  in  this 
scheme  to  a  summation  over  the  points  of  the  fine  mesh.  The 
number  of  terms  in  the  series  being  equal  to  the  number  of 
control  points. 

If  the  pressure  distribution  is  to  be  found  for  a  given 
wing  the  coefficients  of  the  reference  and  additional  lift 
series  are  solved  for  first  and  then  the  pressure  distribution 
can  be  calculated. 

If  an  ideal  wing  camber  is  to  be  defined  from  a  given 
chordwise  pressure  distribution,  the  wing  equation  is  solved 
for  (|rr)p  •  The  terms  in  M*  are  dropped  which  eliminates 
the  leading  edge  singularity,  and  an  equal  number  of  terms  is 
added  on  to  the  end  of  the  series.  In  other  words,  Bos  is 
set  equal  to  zero.  The  induced  drag  of  the  wing  can  be 
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INTEGRATION  POINT 


Fig.  8  Calculation  Meshes 
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reduced  to  an  absolute  minimum  if  the  terms  in  sin  s  e> 
are  eliminated  retaining  only  the  term  in  sin  .  This 
defines  an  elliptical  spanwise  load  distribution.  These  are 
desirable  features  to  incorporate  into  the  design  of  an 
ideal  wing,  and  it  also  simplifies  the  numerical  solution. 

B.  MATRIX  FORMAT 

By  analogy  with  equations  (7.7)  and  (7.8)  the  reference 
and  additional  lift  distributions  may  be  written  in  matrix 
format . 


■l  -  l>]  H 

(8.1) 

J [a] 

,(8.2) 

The  matrix  format  of  the  basic  wing  equations  (6.27)  and  (6.28) 
become 


("If  ~*rj  -  (8-3) 

=  (8.4) 

The  matricies  TqI  and  in  the  above  equations  are  known 

T  by  T  matricies  whose  elements  depend  only  on  wing  planform 
and  aspect  ratio. 

Solving  equations  (8.3)  and  (8.4)  for  the  initially  unknown 
constants  B  and  A  gives  by  matrix  inversion 

—  j 


(8.5) 

(8.6) 
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Equation  (8.6)  shows  that  the  A  constants  depend  only  on 
wing  planform  and  aspect  ratio  while  the  B  constants  in 
equation  (8.5)  also  depend  on  the  wing  slope  function. 

The  points  used  in  the  Q  matrix  need  not  be  the  same  as 
the  control  points  used  in  the  R  matrix,  but  may  be  any 
point  on  the  wing.  In  this  analysis  the  points  for  the  Q 
matrix  were  picked  as  six  stations  across  the  semispan  with 
eighteen  chordwise  stations  at  each  spanwise  station.  This 
permits  a  good  representation  of  the  functions  T  and  (^“)  ^ _ 

If  the  reference  pressure  distribution  over  the  wing  is 
specified,  coefficients  of  the  lift  distribution  series  can 
be  found,  by  equation  (8„7)o  Then  the  wing  equation  (8.8) 
can  be  solved  by  direct  numerical  integration. 


(8.7) 


b*  *r  “  it^ 


TT  IT 


Z  Sr  s'M®  (8.8) 


o  o 


7L 

TT 


+  G,  z  'E  Sr-  Hr-Co^©j  +  T-£g<3B  \  co<^<a  cl© 
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IV  RESULTS 


A.  SOLVING  FOR  THE  WING  PRESSURE  DISTRIBUTION 

The  numerical  results  of  the  pressure  distribution  solution 
were  not  satisfactory.  The  pressure  distributions  and  circu¬ 
lation  functions  calculated  by  the  computer  program  were  not 
smooth  or  in  accordance  with  [Ref.  l]  and  [Ref.  10]  .It  is 
believed  this  is  a  result  of  difficulties  with  the  first 
chordwise  term,  Bos.  The  term  by  term  downwash  is  plotted  in 
Fig.  9  ,  at  a  spanwise  station  near  midspan  for  a  wing  of 
elliptical  planform  and  an  aspect  ratio  of  20  using  5  chord- 
wise  stations  and  10  stations  over  the  span.  The  integration 

mesh  was  four  times  finer  than  the  control  point  mesh.  The 
first  term  should  uniform  downwash.  the  rest  of  the  terms 

are  of  the  general  nature  expected. 
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B.  SOLVING  FOR  THE  WING  CAMBER 

The  numerical  results  of  calculating  the  camber  of  an  ideal 
wing  given  the  desired  chordwise  pressure  distribution  were 
satisfactory.  Fig.  10  compares  the  results  of  — xr  for  a 

wing  of  aspect  ratio  20  and  rectangular  planform  using  a  15 
term  lift  distribution  series  at  5  chordwise  points  near  mid¬ 
span,  with  the  two  dimensional  results  of  [Ref.  10.]  It  is 
felt  the  differences  at  the  leading  edge  can  be  resolved  by 
using  more  terms  in  the  lift  distributions  series  and  more 
integration  points.  The  depicted  pressure  distribution  was 
imposed  at  all  spanwise  stations.  The  results  are  normalized 
to  a  lift  coefficient  of  one.  Table  1  shows  the  results  at 
various  spanwise  stations. 

Table  1 

AR=20  15  term  lift  distribution  series  rectangular  planform 


7 

0.99 

0.89 

0.71 

0.45 

0.16 

i 

^>21 

\z. 

"oZ. 

0.025 

0.0921 

0.1400 

0.1669 

0.1782 

0.1817 

0.206 

0.0292 

0.0561 

0.0789 

0.0883 

0.0913 

0.500 

-  0 . 0443 

-0.0781 

-0.0965 

-0.1038 

-0.1611 

0.794 

-0.0967 

-0.1552 

-0.1951 

-0.2152 

-0.2297 

0.976 

-0.0998 

-0.1456 

-0.1 884 

-0.2142 

-0.2256 

Fig.  10  Camber  Line 
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V.  CONCLUSIONS 


The  numerical  implementation  of  the  theory  needs 
additional  work  to  achieve  the  correct  computer  solutions 
for  the  pressure  distribution  solution.  The  areas  of 
greatest  interest  are  the  assumed  series  solutions  for 
the  additional  and  reference  lift,  especially  the  first 
chordwise  term,  Bos.  The  camber  line  solution  is  consistent 
with  two  dimensional  theory  near  midspan  of  a  high  aspect 
ratio  wing  and  allows  the  camber  line  to  be  determined  at 
various  spanwise  stations  for  any  aspect  ratio.  It  is 
felt  that  this  theory  has  merit  and  deserves  further  study 
to  correct  the  numerical  solution  of  the  pressure 
distribution . 
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APPENDIX  A.  DERIVATION  OF  SUBROUTINES 


A.  SUBROUTINES  FOR  THE  PRESSURE  DISTRIBUTION  SOLUTION 

QMAT  calcualtes  the  Q  matrix  of  equations  (8.1)  and  (8.2) 


s 


1 5>-fe  Sl¬ 


it JV 


to 


Hhx  = 


Z  d  Hr 


d  qS 


Si  =  S/aJ  S  e- 


(7.7) 


(A .  1  ) 


(A. 2) 


(A. 3) 


£>  (  H  s)  =  ±3.  .  |-|  h  1 »  Si. 

c 


(A. 4) 


RMAT  calcualtes  the  R  matrix  of  equations  (8.3)  and  (8.4) 


TT  *' 


(IfL  — *C  = 


TTA^ 

T 

_L_r 

rr 

A 


J  J, 
o  o 


c&,tf  +G,i  IS  w* d 


s 


=  £IB  r*^(Hr)s.MSe 

=  £Z  5rs  H 

dP 


\  n*  ^  s 

Q1-  -  V  r  I-U  S  oot  so 


*  s 

-  3L  ^  (2.  Bos  +  B  <  s  )  S  cot  S 


(6.27) 


(A.  5) 


(A. 6) 


.(A. 7) 


d  « 
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=  i£  Brsj 

“  TT  ii 

TV 

■[Gi 

f-  Q  2.  S  C.oS>  s^j 

ij " 

TT 

tBs]  2  8^  (A  .  8) 

0 

HH1=  ^  (Hr\ 

(A. 9) 

.Si  -  S© 

(A. 10) 

hi  hi  2.  —  f-j  r 

(A. 11) 

S  Z  =  S  COS  S  & 

(A. 12) 

[g.-Hhi-sl  +  c,z-  H(E-d 
0  ^  Ti  AC  -1 

■+■  G>Z  •  S2.  •  J  e>  (A. 13) 

C5 

B.  SUBROUTINES  FOR  THE  CAMBER  LINE  SOLUTION 

QMAT  calculates  the  Q  matrix  of  equation  (8.1)  which  has 
been  simplified  with  a  single  term  in  theta,  an  elliptical 
spanwise  lift  distribution,  and  the  first  chordwise  term  is 
dropped . 


= 

4-^ 

<L 

IS 

d  Hr  1 
dd  -J 

<>l  si  & 

(A. 14) 

l-l  HI  = 

Z 

dHr 

(A. 15) 

d 
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(A. 16) 


<p(  r)  =  .  H  H  1  ♦  Si 

c. 

CAMBER  is  a  numerical  integration  of  equation  (6.27) 
simplified  for  an  ideal  wing. 


bx  TT 


2  I  [>  L 


s/^e 


o  o 


^  1  r- 

+  £  Sr  Hr-  Cos.©!  +  A-  Cote  de 


(A 


.17) 
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